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Abstract 

We investigate hadron attenuation in deep-inelastic lepton scattering off complex nuclei in the 
kinematic regime of the HERMES experiment. Our transport theoretical simulations reveal strong 
prehadronic final state interactions of the reaction products with the surrounding nuclear medium 
early after the initial 7*iV interaction has taken place. In this work we compare our model results 
with the measured hadron multiplicity ratios for a 84 Kr target at HERMES and provide an extended 
discussion of the double-hadron attenuation recently observed at HERMES. 



1 Introduction 

Hadron production in high energy lepton scattering off nuclear targets provides an ideal tool to study 
the physics of hadronization 1 . In elementary electron-proton reactions nothing can be learned about 
the space-time picture of hadron formation, since the reaction products hadronize long before they reach 
the detector. In nuclear reactions the nucleus serves as a kind of micro-detector which is located directly 
behind the primary 7* N interaction point. The attenuation of produced hadrons in the nucleus compared 
to a proton or deuterium target then yields combined information on the length scales of hadron formation 
and the strength of the prehadronic final state interactions (FSI). 

In the recent past the HERMES collaboration jSJEl has carried out a detailed study of hadron produc- 
tion in deep inelastic lepton scattering on various gas targets using the 27.6 GeV positron (electron) beam 
at DESY. The observed attenuation of produced hadrons has led to two complementary interpretations: 
In Refs. 0] it is assumed that hadronization occurs far outside the nucleus and that the attenuation of 
high energy hadrons is solely caused by a partonic energy loss. Here the struck quark undergoes mul- 
tiple scattering while propagating through the nucleus giving rise to induced gluon radiation. On the 
other hand Refs. [SJEJ[7||H1 are based on the idea that the struck quark forms a color neutral prehadron 
early after the initial "f*N interaction and that the prehadronic FSI modify the multiplicity spectra of 
hadrons. The authors of Refs. treat the FSI purely absorptive and restrict their considerations to 

direct interactions of the virtual photon with valence quarks inside the nucleus. In Refs. pillUj we have 
developed a method to incorporate quantum mechanical coherence length effects observed in high-energy 
photon-nucleus reactions into our semi-classical BUU transport model. This enables us to account for 
direct and resolved photon-nucleus interactions and to provide a probabilistic coupled-channel treatment 
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of the (pre-)hadronic FSI. Furthermore, our Monte Carlo approach allows for realistic event-by-event 
simulations that can either be corrected for detector efficiencies or which can directly be used as input 
for detector simulations. In Refs. H3E1 we have applied our model to hadron attenuation at HERMES 
and EMC energies and in Ref. we also made first predictions for the kinematic regime of the current 
JLab experiments. 

Besides making an essential contribution to our understanding of hadronization, electroproduction 
off complex nuclei also provides useful input for the interpretation of observables in ultra-relativistic 
heavy ion collisions. In the latter case the system passes through various energy and density phases 
and the collision geometry is much less under control than for lepton-nucleus scattering where the target 
remains more or less in its ground state throughout the reaction. The hadron formation lengths extracted 
from electron-nucleus reactions help to clarify to what extend the recently observed quenching of high 
transverse momentum hadrons in ultra-relativistic heavy ion collisions at RHIC is due to partonic energy 
loss in a quark gluon plasma EH . 

2 Model 

In our model we split the lepton-nucleus reaction into two parts: In the first step the exchanged virtual 
photon interacts with a bound nucleon inside the nucleus, thereby producing a prehadronic final state 
which - in step two - is propagated in our BUU transport model. 

In the initial 7*iV interaction we account for nuclear effects such as binding energies and Fermi motion 
of the bound nucleons as well as Pauli blocking for outgoing fermions. Furthermore, we calculate the 
wave function modifications of the resolved photon components, i.e. the vector meson and perturbative qq- 
fluctuations, on the path through the nucleus up to the primary production point. As a result the nucleons 
on the front side of the nucleus shadow the downstream nucleons from resolved photon interactions at 
large coherence lengths The situation is depicted in Fig. ^ for a 84 Kr nucleus and a high-energy 

virtual photon with momentum directed along the positive z-axis and a p° coherence length of the order 
of the nuclear radius. 

The product of the high-energy j*N interaction is determined by the event generators Pythia JJ] 
and Fritiof JH] which allows us i) to account for direct interactions of the bare virtual photon with sea 
and valence quarks inside the nucleon and ii) also for the resolved photon interactions. In the latter case 
the hadronic fluctuations of the photon can participate in a hard or a diffractivc scattering process with 
the nucleon. All scatterings lead to the excitations of hadronic strings which subsequently decay according 
to the Lund fragmentation scheme |16| into color neutral prehadrons. In this work we assume that these 
prehadrons are produced directly after the j*N interaction. This is a simplifying assumption that is 
generally used in transport models 117| . The prehadrons then need a formation time to build up their 
hadronic wave function and to interact like physical hadrons. For simplicity we take the formation time to 
be a universal constant tj — 0.5 fm/c in the rest frame of each hadron. This corresponds approximately to 
the time that the hadron constituents need to travel a distance of a typical hadronic radius. Furthermore, 
it is in agreement with formation times that have been extracted from particle production in pA collisions 
at AGS energies within a similar transport model 18 . During their formation time the prehadrons are 
assumed to interact with reduced cross sections which are determined by the number n org of leading 
(anti-)quarks that were initially present either in the target nucleon or the hadronic fluctuation of the 
photon: 

* _ n Q'g 
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Figure 1: Profile function for shadowing. Shown is the probability distribution for the interaction of an 
incoming photon (from left) with given virtuality Q 2 and energy v with nucleons in a 84 Kr nucleus. 
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Consequently, the hadrons that contain leading quarks, i.e. the beam and target remnants, interact 
with a (reduced) cross section right after the j*N interaction while the prehadrons that solely contain 
(anti)quarks produced from the vacuum in the string fragmentation will start to interact after their 
formation time tf — jTf. 

The propagation of the final state (pre) hadrons through the nucleus is described by our BUU transport 
model. The model is based on a set of generalized transport equations which determines the time evolution 
of the spectral phase-space distributions Fh of the (pre)hadrons. The distributions Fh(r,p,/j>,t) give at 
each moment of time t and for each particle class h the probability to find a particle of that class with a 
(possibly off-shell) mass p, and momentum p at position r. Its time-development is given by the equation 



d dH h d dH h 8 



dt dp dr dr dp 



F h = G h A h - L h F h (2) 



where Hh denotes the relativistic energy of the (pre)hadron h and contains a momentum- and density- 
dependent mean-field potential in case of baryons. The terms on the lhs of (01 are the so-called drift terms 
since they describe the independent transport of each hadron class h. The terms on the rhs of @ are the 
collision terms; they describe both elastic and inelastic collisions between the hadrons. Here the term 
inelastic collisions includes those collisions that either lead to particle production or particle absorption. 
The former is described by the gain term GhAh on the rhs in J2J, the latter process (absorption) by the 
loss term LhFh- Note that the gain term is proportional to the spectral function Ah of the particle being 
produced, thus allowing for production of off-shell particles. On the contrary, the loss term is proportional 
to the spectral phase space distribution itself: the more particles there are the more can be absorbed. 
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The terms Gh and Lh on the rhs give the actual strength of the gain and loss terms, respectively. They 
have the form of Born-approximation collision integrals and take the Pauli-principle into account. The 
free collision rates themselves are taken from experiment or are calculated EH . 



3 Results 

In Ref. [7] we have presented a detailed comparison of our model results and experimental data from the 
HERMES experiment. In this work we review our major findings. 

In our simulations we apply the same kinematic cuts as in the HERMES experiment and in addition 
account for the geometrical acceptance of the HERMES detector. The first observable that we discuss is 
the so called multiplicity ratio 

R M (z h ,v,p T )- Nh{Zh ^ Th , (3) 
n7(7) Id 

where Nh is the yield of semi- inclusive hadrons in a given (zh, i/,Py)-bin and N e the yield of inclusive 
deep inelastic scattering leptons in the same i/-bin. Like in experiment we define the hadron fractional 
energy as Zh = Eh/v and denote the transverse hadron momentum component with respect to the virtual 
photon momentum as px- For the deuterium target, i.e. the denominator of Eq. we simply use the 
isospin averaged results of a proton and a neutron target. Thus in the case of deuterium we neglect 
the FSI of the produced hadrons and also the effect of shadowing and Fermi motion. The latter will be 
accounted for in an upcoming work. 

Fig. [3 shows the multiplicity ratio for the production of pions, (anti-)kaons and (anti-)protons 
on a 84 Kr target. The solid line represents the result of our coupled-channel BUU calculation using an 
universal formation propertime t/ = 0.5 fm/c and our simple constituent quark concept As we have 
discussed in Ref. the increase of with photon energy v is not solely caused by an increase of 
time dilatation but may in parts be due to the geometrical acceptance of the HERMES detector which is 
restricted to the forward region. The increase of hadron attenuation with rising fractional hadron energy 
Zh arises from the fact that a high Zh (pre) hadron most likely contains a leading quark and therefore 
might get absorbed early after the 7* A interaction. Furthermore, the elastic and inelastic FSI shuffle 
strength from higher to lower (fractional) energies as can be seen by comparing to a calculation without 
coupled-channel treatment (dashed line). In the latter case we simply remove all particles that have 
undergone a collision from the outgoing channel. As expected the discrepancy is largest at small Zh- One 
also sees that the coupled-channel treatment leads to an increasing yield of low energy protons that are 
knocked out of the nucleus in the FSI. However, the kinematic cuts reduce this increase in such a way 
that we still overestimate the attenuation of protons. 

From Fig. [3 one sees that the absorption strength is about the same for kaons and anti-kaons. This is 
due to the following effects: Because of reactions like KN — ► Air . . . the absorption cross section of K~ is 
larger than that of K + . This is partly compensated by a reduction of the prehadronic interactions. This 
is due to a smaller probability to form a leading K~(us) in the primary 7* A interaction, since now the 
photon has either to interact directly with a sea quark inside the nucleon or undergo a resolved interaction. 
Furthermore, there exists the possibility to create K + {us) through reactions like irN — > A K . . . in the 
FSI. 

The HERMES collaboration has also measured the pr-dependence of the multiplicity ratio © . The 
result is shown in Fig. [3] for a 84 Kr target in comparison to our BUU calculation indicated by the solid 
line. Obviously, the experimental multiplicity ratio shows a strong increase at transverse momenta above 
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Figure 2: Multiplicity ratios of tt*' , K^, p and p for a 84 Kr nucleus (when using a 27.6 GeV positron 
beam at HERMES) as a function of the hadron energy fraction Zh = E-^jv and the photon energy v. 
The solid line represents the result of a simulation, where we use the constituent quark concept for 
the prehadronic cross sections and a formation time t/ = 0.5 fm/c. The dotted line indicates the result 
of a simulation where the FSI have been treated purely absorptive. The data are taken from Ref. [Hj- 

1 GeV which can be seen as the Cronin effect in eA collisions and which is not reproduced by our 
(pre)hadronic FSI. In addition to our standard approach - which uses the vacuum cross sections as input 
for determining the outcome of the rescatterings in the FSI - we have also performed a calculation where 
we assumed all elastic scatterings in the FSI to be isotropic (dashed line in Fig. |3J). The latter leads 
to an increase of the theoretical multiplicity ratio at high transverse momenta, however it overestimates 
the total strength of the attenuation! The Cronin effect in eA collisions therefore should have a partonic 
origin which might cither be due to partonic rescattering in the FSI, an increase of the intrinsic transverse 
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Figure 3: The multiplicity ratio of charged hadrons for 84 Kr (at HERMES) as a function of the transverse 
momentum squared p^. In the simulation we use the constituent quark concept Q for the prehadronic 
cross sections and a formation time Tf = 0.5 fm/c. In the simulation indicated by the dashed line we 
additionally assume that all clastic scatterings are isotropic in their center-of-mass frame. The data are 
taken from Ref. py. 



momentum {kr) of the resolved photon components on their way through the nucleus - an initial state 
effect - or an increase of (fcy) in the bound target nucleons. 

Finally, we discuss the double-hadron attenuation that has recently been measured at HERMES 20 . 
In each event only the two (charged or neutral) hadrons with the highest energies are considered. In 
the following we denote the hadron with the highest Zh as the leading hadron and the other one as the 
subleading hadron. The experimental observable is the double-hadron attenuation ratio 



R2(z 2 ) = 



N 2 (z 2 ) 


A 


Ni 


N 2 (z 2 ) 


D 


Ni 



(4) 



Here ^2(^2) denotes the number of events where the leading and subleading hadron carry the energy 
fraction Z\ > 0.5 and z-i < z\, respectively, and N\ is the number of events where at least one of them has 
Zh > 0.5. Again we account for all experimental cuts and the geometrical acceptance of the HERMES 
detector in our simulation. 

Fig.0Jshows our model prediction for the double-hadron multiplicity ratio (0} for 14 N, 84 Kr and 131 Xe. 
To exclude contributions from p° decay into 7r + 7r~ the charge combinations 'H — ' and '— + ' have been 
excluded both in experiment and in the simulation. The solid line shows the result of a full coupled- 
channel calculation using the constituent quark concept Q and the formation time 77 = 0.5 fm/c. The 
shape of the spectrum is similar to that of R^(zh) of the single hadron multiplicity ratio shown in Fig. [21 
As can be seen from Fig. QJ our calculations predict about the same double-hadron attenuation ratio for 
131 Xe and 84 Kr. As we explain below the attenuation of leading and double hadrons increases in the same 
way when going from the krypton to the xenon target. Hence, the double-hadron attenuation ratio 
stays roughly the same. Note, that this does not necessarily imply the same single-hadron attenuation for 



(i 



84 Kr and 131 Xe. As we have shown in Refs. [HIE] the single-hadron attenuation @ in the 131 Xe nucleus 
is larger than for 84 Kr. 

To understand the behavior of the double-hadron multiplicity ratio qualitatively we consider a 
simple absorption model (cf. [5]). We assume that the leading prehadrons that contribute to the ratio 
Ri = Ni\a/Ni\d start to interact with a cross section a at distance l p behind the 7*iV interaction point 
(b,z), i.e. 

oo oo 

Ri = — J d 2 b J dz pA(b, z) exp [ — a J dz' pA(b, z')\ , (5) 

— oo z+l p 

where A denotes the mass number of the target nucleus and pa its density. If one assumes the production 
lengths l\ and 1% for the leading and subleading prehadrons, one gets for the double-hadron multiplicity 
ratio 

oo oo oo 

R-i = ^ J d 2 b J dz pA(b 1 z) exp [ — a J dz' p A {b 1 z')] exp [ — a J dz' p A (b, z')] . (6) 

-oo z+h z+l 2 

Even if all production lengths are the same (l p = li = 1% = V), R2 will differ from i?i: 

^ = / d 2 b dz p A (b, z) exp [ - 2a dz'p A {b, z')} ^ 
J d 2 b dz p A (b, z) exp [ - a dz'p A (b, z')] 

The numerator and the denominator in Eq. Q differ by the absorption cross section in the exponent 
which for the case of double-hadrons is twice the cross section for the single hadron absorption. Hence, 
i?2 is smaller than one for pure absorption and depends only very weakly on the mass number A of the 
target nucleus, since the A-dependence of the numerator and the denominator are about the same. The 
dashed line in Fig.0|shows the result of our BUU model for a purely absorptive treatment of the FSI. The 
increase of R2 (2^) with decreasing z-i is obviously due to the coupled-channel effects: For the interpretation 
we discard for a moment the constant factors Ni in Eq. |0J and the factors N e in Eq. © which have 
no influence on the shape of the Zh dependence. The only difference between Nh{zh)\A/Nh{zh)\D and 
N2{z2)\a/N2{z2)\d then is that one restricts the detected hadron to the subleading particle in the latter 
case. If the subleading particle (with energy fraction Z2) of the initial 7*^ reaction interacts with the 
nuclear environment, it will produce a bunch of low-energy particles. The new subleading hadron in the 
event then has a energy fraction z' 2 < Z2- As for the usual charged hadron multiplicity spectrum the 
coupled-channel FSI shuffle strength from the high Zh part to the low Zh part of the spectrum. This 
is not the case for purely absorptive FSI (dashed line in Fig. QJ. As one can see, our coupled-channel 
calculations (solid line) - using the constituent quark concept Q and the formation time Tf = 0.5 fm/c 
- are again in quantitative agreement with the experimental data apart from the last data point in the 
84 Kr data, which indicates a multiplicity ratio i?2(-^2 = 0.5) « 1. This behavior cannot be explained 
within our model. 



4 Summary &; Outlook 

We have presented a model which combines a quantum mechanical treatment of coherence length effects 
in high-energy (virtual) photon-nucleus interactions with a probabilistic coupled- channel description of 
hadronic FSI. This allows us to model hadron electro- and photoproduction off complex nuclei at multi 
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Figure 4: Double-hadron attenuation ratio R2 for a N, 8 Kr and 131 Xe target as a function of the energy 
fraction Z2 of the subleading hadron. In the simulation (solid line) we use the constituent quark concept 
P| for the prehadronic cross sections and a formation time 77 = 0.5 fm/c. To exclude contributions from 
p decay into tt + w~ the charge combinations 'H — ' and ' — h' have been excluded both in experiment and 
in the simulation. The dashed line shows a calculation with a purely absorptive treatment of the FSI. 
The preliminary HERMES data are taken from Ref. [20) . 



GeV energies in a way that goes far beyond simple absorption models. Furthermore, the model can be 
used with basically the same parameters and physics assumption to describe a broad range of different 
nuclear reactions such as pion and proton induced reactions as well as heavy-ion-collisions. 

In this work we have compared our model results to experimental data on single- and double-hadron 
multiplicity ratios in deep inelastic positron scattering off nuclear targets at HERMES. We achieve a 
good description of most experimental data by assuming early interactions of color neutral prehadrons 
that are produced in the initial 7* TV interaction and account for particle absorption and production in 
the FSI. However, the model underestimates the increase of high transverse momentum hadrons in eA 
reactions which may indicate additional partonic effects in the FSI. 

We are currently optimizing our model in the HERMES kinematic regime by implementing radia- 
tive corrections and Pythia parameters that have been adjusted by the HERMES collaboration \'21\ . 
Furthermore, we will extend our studies to eA scattering at JLab (221 7 EMC |231 an d EIC [21] energies 
in the future. Since our model works on an event-by-event basis its output can later on be used for 
various detector studies. In addition we will improve our space-time description of hadronization. In 
Ref. we have developed a method to extract the complete four-dimensional space-time coordinates 
of the prehadron production and hadron formation points from the Jetset routines of Pythia on an 
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event-by-event basis. This information will be implemented into our BUU transport model. Furthermore, 
one can think of considering an additional partonic energy loss prior to color neutralization in our model. 
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